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Abstract. We study the cooling of isolated neutron stars. The main cooling regulators are introduced: EoS, 
thermal transport, heat capacity, neutrino and photon emissivity, superfluid nucleon gaps. Neutrino emissivity 
includes main processes. A strong impact of medium effects on the cooling rates is demonstrated. With taking 
into account of medium effects in reaction rates and in nucleon superfluid gaps modern experimental data can be 
well explained. 

Key words, dense baryon matter, neutron stars, medium effects, pion softening, nucleon gaps, heat transport 



1. Introduction 

The Einstein Observatory was the first that started 
the experimental study of surface temperatures of iso- 
lated neutron stars (NS). Upper limits for some sources 
have been found. Then ROSAT offered first detections 
of surface temperatures. Next X-ray data came from 
Chandra and XMM/Newton. Appropriate references to 
the modern data can be found in recent works by 
ITsuruta et al. 2002L ITsuruta 20041 IKaminker et al. 20011 
lYakovlev at al. 2003al devoted to the analysis of the new 
data. More upper limits and detections are expected from 
satellites planned to be sent in the nearest future. In gen- 
eral, the data can be separated in three groups. Some data 
show very "slow cooling" of objects, other demonstrate an 
"intermediate cooling" and some show very "rapid cool- 
ing". Now we are at the position to carefully compare the 
data with existing cooling calculations. 

2. Existing NS cooling scenarios 

Let us briefly point out some important achiements on the 
way to the present understanding of the NS cooling prob- 
lem. Theoretical study of the NS cooling has been started 
long ago in pioneering works of Tsuru ta fc Cameron 1965 



and lBahcall fc Wolf 19651 It has been argued that the one- 
nucleon so called direct Urea (DU) process, as n ^ peD, 
is forbidden up to sufficiently high density and the main 
role plays the two-nucleon so called modified Urea (MU) 
process, like nn npev and np ppev. As the result of 
many works the so called "Standard scenario" of cooling 
emerged. It includes the neutrino cooling stage t ^ 10^ yr 
and the photon cooling era, t Z 10^ yr. The MU pro- 
cess and the nucleon- nucleon bremsstrahlung (NB), as 
nn —>■ nnvv and np — » npi/v, were carefully recalcu- 
lated bv lFriman <k Maxwell 19791 using the free one-pion 
exchange. The emissivity of the most efficient channel of 
the MU process nn — > npeu is given by 
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Tg ~ T/IO^K, rig is the electron density and uq ~ 
0.16 fm~^ is the ordinary nuclear density at saturation, 
('^Mu("))^ — (?Ti*i/"iAr)"^(r7i*/TOAr), m* and m* are effec- 
tive neutron and proton masses, m^r is the nucleon mass in 
vacuum, n is the nucleon density, {m*j^j-^{n))^{ne/nQy^'^ ^ 
0.1 for n ~ no. 

After [Migdal 195*91 11 became clear that NS in the late 
cooling stage for T < Ten and T < T^p are neutron and 
proton superfluids. The nucleon superfluidity was incor- 
porated in the "Standard scenario " bv lMaxwell 1971 who 
used relevant combinations of the suppression factors 



~ exp(-Aii/T), T < Tc, 
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for the emissivity of MU nn 
■1 



n-pev and NB np 

2 



n-pvy processes, for NB nn — > nnvD, S^^^ for MU np — > 
ppev, for T < Ten and T < Tcp, A,;i is the gap, i = n or p. 

Typically the NB emissivity is by an order of magni- 
tude smaller than that for MU. The pre-factor ~ 10^° 
10^^ and the temperature dependence Tg are typical for 
the phase space volume of two-nucleon processes (without 
inclusion of medium effects, see below). 

First cooling calculations were based on the as- 
sumption of an isothermal core. The relevance of the 
thermal conductivity for the first 10^ 10'^ yr was 
demonstrated by INomoto fc Tsuruta 19811 A discus- 
sion of the achievements of early works can be found 
in the review by Tsuruta 1979. and in the book by 
[Shapiro fc S.A. Teukolsky 1983| 

The "Standard scenario" allows to explain the slow 
cooling but fails to cook the intermediate cooling and 
fast cooling of some neutron stars. To explain the lat- 
ter different efficient direct Urca-like processes have been 
suggested: the proper DU process, reinvestigated by 
ILattimer et al. 19911 and the DU process going on hy- 
perons (HDU), see the same work; pion condensation 
(PU) processes by IMaxwell et al. 19771 kaon condensa- 
tion (KU) processes bv lBrown et al. 19881 iTatsumi 19881 
similar to PU; and quark DU processes (QDU) by 
P^wamotp 19821 All these processes show up for baryon 
densities larger than the corresponding critical ones n > 
nj?u, nPU, nf"^, n™^^ „QDU _ ^ according to dif- 
ferent model dependent estimates. Roughly 
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that this important result was overlooked for many years. 
[Voskresensky fc Senatorov 1987, rediscovered the nPBF 
process n nvv and introduced the pPBF process 
p — > pvD with correct asymptotic behavior of the emis- 
sivity 



where (toqu)" 
and the behavior ^ Tg are typical for all mentioned one- 
nucleon processes. Thus the "Standard scenario + exotics" 
has been considered as a scenario with minimum exotics: 
some stars have n < nc and cool slowly whereas some have 
n > ric in a part of their interiors and cool very fast. 

However one type of the non-exotic processes was com- 
pletely forgotten already in the "Standard scenario" with- 
out any justification for that. These are the so called 
neutron pair breaking and formation (nPBF) and proton 
pair breaking and formation (pPBF) processes, as they 
were named bv lSchaab et al. 19971 The emissivity of the 
nPBF process was first calculated by flowers ct al. 1976| 
for the case of the 15*0 neutron pairing. However, their 
asymptotic expression for the emissivity e[nPBF] ^ 
10^0 T^ exp(-2A„„/r) for T < A™, as follows from 
expression (lb) of their work and from their rough asymp- 
totic estimate of the integral (see below (13b)), shows nei- 
ther the large one nucleon phase space factor (~ 10^^) 
nor the appropriate temperature behaviour (note that 
the full analytic expression of this work is quite cor- 
rect). The numerical underestimation of the rate by an 
order of magnitude and the very rough asymptotic ex- 
pression used by IFlowers et al. 19761 became the reason 
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for T < A,, 
mentum. 



ppi{n) is the neutron/proton Fermi mo- 
i = "i^pgp = m*/mN. Eq. shows 

a huge one-nucleon phase space factor and a quite 
moderate T dependence of the pre-factor (note that 
the value of the pre-factor in case of 15*0 neutron 
pairing evaluated by [Voskresensky fc Senatorov 1987| 
contains a mistake that however yields only a non- 
essential correction factor of order of one to their final 
result; a larger uncertainty comes from not too well 
known medium modified vertices, see discussion by 
[Voskresensky 200T] and corrected expressions there). 
Medium modifications of vertices (cf. Ward-Takahashi 
identity) were incorporated in all processes. This is very 
important point since medium modified vertices yield 
an enhancement factor up to 10^ for pPBF compared 
with the result one would get using vacuum vertices. 
This factor ^ 10^ arises since the process may occur 
through nn^^ and ee~^ correlations, —1 symbolizes the 
particle hole, with subsequent production of from the 
n and e rather than from a strongly suppressed vacuum 
vertex of pi^i', see Voskresensky fc Sena torov 1987', 



Voskresensky et al. 1998[ iLeins on 2000j 

^osl^resens^^^00l | The efficiency of the PFB rates 
was related by I Voskresensky fc Senatorov 1987[ to the 
value of the pairing gap and it was indicated that these 
processes can play an important role in the cooling 
scenario. Note that with taking into account medium 
effects the estimate (0)) is roughly valid for all the 
cases for ISq neutron and proton pairings as well as 
for the 3P2 neutron pairing, except the specific case, 
when the projection of the pair momentum onto the 
quantization axis is \mj\ = 2. In the latter case, where 
the gap has zero's in some points of the Fermi surface, 
the emissivity would not be exponentially suppressed, 
cf. [Voskresensky fc Senatorov 1987[ However conditions, 
when such a possibility might be realized, are not clear. 

The inclusion of the PBF processes in the cooling code 
has first been accomplished by ISchaab et al. 19971 (where 
also for the first time the reference to the pioneering work 
of Flowe rs et al. 1976i appeared) with the result that vari- 
ation of gaps allows to develop a "scenario for intermedi- 
ate cooling" covering the region of both "slow and inter- 
mediate cooling". Some subsequent works rediscovered the 
pPBF process with vacuum vertices thus ignoring medium 
effects but including relativistic corrections to the vacuum 
vertex. The latter corrections yield an enhancement of or- 
der of 10 rather than the mentioned 10^ and are therefore 
of minor importance compared to the larger contribution 
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of medium effects. The pPBF emissivity evaluated with 
vacuum vertex (with relativistic corrections) was imple- 
mented into their cooling code bv lYakovlev at al. 2003al 
Their result for this process should be still enhanced at 
least by a factor > 10 to account for medium effects. 

Medium effects essentially modify not only 
the pPBF emissivity but they correct also the 
contributions of all other processes. They have 
first been included into emissivities of differ- 
ent processes by IVoskresensky & Senatorov 1984| 



Voskresensky & Senatorov 1986 


see also 


Voskresensky & Senatorov 1987 


Migdal et al. 1990| 


and a 


more recent review 


Voskresensky 2001 


It was shown that the main contri- 



bution to the MU process actually comes from the pion 
channel of the reaction nn — > npev, where ev are radiated 
from the intermediate pion and the NN~^, exchanging 
nucleons, rather than from the nucleon of the leg of the 
reaction. Moreover, due to the so called pion softening 
(medium modification of the pion propagator) the matrix 
elements of the MU process are further enhanced with 
the increase of the density towards the pion condensation 
critical point. Thus the corrected MU process has been 
called "the medium modified Urea" (MMU) and corrected 
NB process was called MNB. Roughly, the emissivity 
(|T]l acquires then a factor (mainly due to the pion decay 
channel of MMU) 



e^[MMU] 



10^ (n/no)i°/^- 



e.[MU] 
where the pre-factor 



(5) 



(n/rio)^"^'^ arises from the phase 
space volume, T{n) = 1/[1 -(- C{n/nQf/% C ~ 1.4-^-1.6, 
is the proper nucleon- nucleon correlation factor, which 
appears in the vertices due to the strong interaction, 
being expressed through the Landau-Migdal parameters 
(in another words r(n) recovers the Ward-Takahashi 
identity). Note that the above introduced value F is an 
averaged quantity. Actually correlation factors depend 
on the energy-momentum transfer, being different for 
vertices connected to the weak coupling (F^j-g is rather 
close to unity) and for vertices related to the pure strong 
coupling (Fs is slightly less than above introduced factor 
F, F^ = F^_gFg). uj*{n) is the so called effective pion gap. 
The latter has the meaning of the effective pion mass at 
finite momentum transfer in the given reaction, k ~ ppn- 
The quantity w*^(n) replaces the value {m^ + pj^^) in the 
case of the free pion propagator used for the calculation 
of the MU process bv lFriman fc Maxwell 19791 Thus the 
squared amplitude of the nucleon-nucleon interaction 
contains ^ F^/w*^ instead of ~ l/lw^ + P'pn] fo'^ 
free one pion exchange used bv lFriman fc Maxwell 19791 
Here we suppressed a smaller contribution of a lo- 
cal Landau-Migdal interaction, being corrected by 
NN~'^ loops, cf. [Voskresensky fc Senatorov 1986| 
IVoskresensky fc Senatorov 1987| [Migdal et al. 1*9901 

The density dependencies of the correlation factor and 
the effective pion gap are presented in Fig. We see 
that vertices are rather strongly suppressed (and this 
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Fig. 1. Nucleon - nucleon correlation factor F and squared 
of the effective pion gap lu* with pion condensation 
(branches la, 2, 3) and without (la, lb). 



suppression increases with the density) but the soften- 
ing of the pion mode is enhanced (cj*^ < to^) for n -J 
0.5 0.8 uq. Such a behavior is motivated both the- 
oretically and by analysis of nuclear experiments, see 
[Migdal et al. 1990| lEricson fc Weise 19881 The curve la 
shows the behavior of the pion gap for n < n^^ . The 
value n^^ depends on different rather uncertain param- 
eters (we further assume n^^ ~ 3 no) 7 cf. discussion by 
[Migdal et al. 1990[ The curve lb demonstrates the possi- 
bility of a saturation of pion softening and the absence 
of pion condensation for n > n^^ (this possibility could 
be realized, e.g., if Landau-Migdal parameters increased 
with the density). Curves 2, 3 demonstrate the possibil- 
ity of pion condensation for n > n^^ . The continuation 
of the branch la for n > n^^ (the branch 2) demon- 
strates the reconstruction of the pion dispersion relation 
on the ground of the condensate state. Here for simplicity 
we do not distinguish between 7r°, n"^ condensations, see 
[Voskresensky fc Senatorov 1984[ [Migdal et al. 1990[ and 
the so called alternative-layer-structure, including both 
types of condensates, see [Umeda et al, 1994[ In agree- 
ment with a general trend known in condensed matter 
physics fluctuations dominate in the vicinity of the crit- 
ical point of the phase transition and die out far be- 
low and above the critical point (see the curves la, 2). 
The jump from the branch la to 3 is due to the first 
order phase transition to the it condensation, see dis- 
cussion of this point by [Voskresensky fc Mishustin 1982| 
[Migdal et al. 1990[ The branch 3 yields the amplitude of 
the pion condensate mean field for n > n^^ . The ob- 
servation that the pion condensation appears by the first 
order phase transition needs a comment. With the first or- 
der phase transitions in the systems with several charged 
species is associated the possibility of the mixed phase, see 
[Glendenning 1992| The emissivity is increased within the 
mixed phase since efficient DU-like processes due to nu- 
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cleon re-scattering on the new-phase droplets are possible. 
However [Voskresensky et al. 2002| [Maruyama et al. 2003| 
demonstrated that, if exists, the mixed phase is probably 
realized only in a narrow density interval due to the charge 
screening effects. Thereby to simplify the consideration we 
further disregard the possibility of the mixed phase. We 
also disregard the change in the equation of state (EoS) 
assuming that the phase transition is rather weak. 

The works of jVoskres ensky fc Senatorov 1984| 
[Voskresensky fc Senatorov 1986| suggested to explain the 
difference in the surface temperatures of various compact 
objects by the assumption that the objects have different 
masses and thus different density profiles and different 
cooling rates, according to above argumentation. 

The most precise measurement of NS masses comes 
from the binary pulsar system PSRB1913+16, where for 
the more massive a value M = 1.4411 ± 0.0007 M© has 
been obtained [Taylor fc Weisberg 1989| The approximate 
coincidence of this value with the present statistical aver- 
age value for masses of binary pulsars M « 1.35±0.04 Mq 
[Thorsett &: Chakrabarty 1999[ motivated some authors to 
conjecture that all NS masses should be very close 1.4 Mq. 
Recent measurement |Lyne et al, 2004| of the NS mass 
M ~ 1.25Mq in JC737-3C3J does not confirm the above 
conjecture. NS masses can vary in some limits. Thus the 
argumentation, cf. [Voskresensky fc Senatorov 1986[ that 
slow to fast cooling is explained by different masses of 
the corresponding objects seems rather natural. Many 
other in-medium reaction channels have been studied. The 
particular role of medium effects in the NS cooling has 
been reviewed by 'Migd al et al. 1990| and more recently 
by [Voskresensky 2001, The pion softening with increas- 
ing baryon density and the subsequent pion condensa- 
tion are now-days reproduced not only by microscopic 
models of the polarization of the baryonic medium, see 
[Migdal et al. 1990[ and lSuzuki et al. 19991 but also by de- 
tailed variational calculations of lAkmal et al. 19981 yield- 
ing ~ 1.5 -^- 3 no for the tt" and charged n~ con- 
densates. We repeat here the conclusion of the series of 
above mentioned works that only due to the enhancement 
of medium polarization with the baryon density the pion 
condensate may set in and it seems thereby not justified 
to ignore the softening effects for n < n^^, and suddenly 
switch on the condensate for n > n^^ . 

In-medium modifications of different processes includ- 
ing MU, NB, PFB, PU and DU reactions were imple- 
mented in the cooling code by Schaab ct al.J^7 with the 
conclusion that the data existing to that time could be 
explained assuming different masses of the sources. The 
statement was demonstrated on different models for the 
equation of state (EoS) with and without superfluidity ef- 
fects. However, the best fit to the whole set of the data was 
not elaborated. E.g., the slow cooling object PSR 1055- 
52 was not fitted. As we show below it cannot be done 
without a fit of the 15*0 n and p and the 3P2 n gaps, 
their absolute values and density dependencies. An alter- 
native assumption of an internal heating suggested to ex- 
plain a high temperature of this object, although possible. 



seems us rather artificial and we drop it. Inclusion of dif- 
ferent medium effects (see [Voskresensky 2001) ) allows to 
develop a general "Nuclear medium cooling scenario " de- 
scribing all the cooling data and covering the gap between 
"Standard" and "Standard+ exotics" scenarios. The works 
of IBlaschke et al. 20001 and IBlaschke et al. 20011 consid- 
ered possibility of color superconducting pure quark and 
hybrid stars discussing peculiarities of these objects. 

The newly appeared data were studied in recent 
works by Kaminkcr ct al. 2001' "Yak ovlev at al. 2003al 
and .Tsuru ta ct al. 2002 , Tsuruta 20Q4| Medium ef- 
fects in the emissivity were disregarded. The cooling 
scenario based on the inclusion of the DU process 
( "Standard+ exotics DU") has been worked out by 
IKaminker et al. 20011 lYakovlev at al. 2003al Their result 
is that the neutron 3P2 gap should be strongly suppressed 
since otherwise one could not explain the "slow cooling" 
objects. This assumption is supported by some modern 
microscopic calculations of the 3P2 gap including medium 
effects, see .Schulze et al. 1996 Lombardo & Schulze 20001 
ISchwenk & Friman 20081 Contrary, according to 
IKaminker et al. 20011 lYakovlev at al. 2003a] the 15*0 
proton gap should be enhanced in order to smoothen the 
transition from the "slow cooling" to the "rapid cooling". 
Only switching on the DU process for M > 1.358M0 
in their EoS model allows to explain the " intermediate 
cooling" and "rapid cooling" of some objects. The 
transition from " slow cooling" to "rapid cooling" occurs 
in a narrow window of masses. 

Criticism of the Yakovlev et al. DU "Standard+exotics 
DU" scenario is as follows, see also ITsuruta et al. 20021 
ITsurut a 2004 i) The most elaborated model of the vari- 
ational theory EoS of lAkmal et al. 19981 permits the DU 
process only for rather high density n > ~ 5.19 no 
and M > 1.839 AIq, whereas lYakovlev at al. 2003al used 
a model of the EoS, where the DU process starts for 
n > ~ 2.8 no (M > 1.358 Mq). It seems doubtful 
that the DU reaction "might know" that it should start 
namely, when the NS mass is approximately 1.4 Mq, be- 
ing very close to the average value of the NS mass mea- 
sured in NS binaries. The NS mass is governed by the 
full NN interaction whereas the value of the proton frac- 
tion, being responsible for the critical density of the DU 
reaction, is governed by only a part of the interaction re- 
lated to the symmetry energy. In the light of the recent 
observation of the object with M ~ 1.25 Mq the above 
assumption seems even more unlikely. If the DU process 
occurred at M > 1.839 Mq as follows from the EoS 
bv lAkmal et al. 19981 the model of lYakovlev at al. 2003al 
based on the opening of the DU channel would fail. 
Indeed, it would mean that the majority of the isolated 
NS have masses M > 1.839 Mq. The latter sounds quite 
unrealistic, ii) lYakovlev at al. 2003a] introduced by hand 
a large proton gap. However, there are indications by 
ISchulze et al. 19961 that medium effects should suppress 
the proton gap. Moreover, [Takatsuka fc Tamagaki 19971 
demonstrated that proton and neutron superfluidities 
must be very weak for densities when the DU process be- 
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comes operative, iii) As we argued, the pPBF cmissivity 
should be enhanced by an order of magnitude by medium 
effects, what is not taken into account. Moreover, iv) the 
MU and the NB processes are significantly affected by 
medium effects (should be replaced to MMU and MNB) 
what is not incorporated. 

Works of ITsuruta et al. 20021 TTsuruta 200l used the 

"Standard+ exotics PU" scenario (based on the inclusion 
of the PU process for n > n^^ ~ 2.5 uq and Auq respec- 
tively). In general their "slow cooling" curves are below 
those of lYakovlev at al. 2003al Thereby the hottest object 
PSR 1055-52 (if the surface temperature of it is correctly 
extracted from the measurement) is not described, similar 
to the result of .Schaab et al. 19971 The transition from the 

"slow cooling" to the "rapid cooling" is here due to the 
PU process. Again the model assumes that the value of the 
transition density for the non-standard process (here pion 
condensation) is very close to the value of the central den- 
sity of a "magic" star of M ~ 1.4 Mq in order to have all 
NS masses very near the value of 1.4 Mq. The latter seems 
rather unlikely in the light of the recent measurement of 
the NS star mass M ~ 1.25 Mq, see |Lyne et al, 2004| The 
medium effects are disregarded and thus the above points 
iii), iv) remained not corrected. 

The necessity to include in-medium effects into the NS 
cooling problem is a rather obvious issue. It is based on 
the whole experience of condensed matter physics, of the 
physics of the atomic nucleus and it is called for by the 
heavy ion collision experiments, see [Migdal et al. 1990| 
IVoskresensky 2001| llvanov et al. 200 11 111 Their rele- 
vance in the NS cooling scenario was demonstrated by 
ISchaab et al. 19971 Below we present cooling calculations 
of NS based on the "Nuclear medium cooling scenario", 
see IVoskresensky 200T] We remove the above shortcom- 
ings and show the relevance of medium effects on the 
reaction rates and superfluid gaps by demonstrating the 
possibility to fit the whole set of cooling data available by 
today. 



3. EoS and structure of NS interior, crust, surface 



the original Argonne V18 + 5v + UIX* model yielding 



5.19 no (Mpu ~ 1.839 Mq). 



3.1. NS interior 



We will exploit the EoS of lAkmal et al. 19981 (specifically 
the Argonne V18 + 5v + UIX* model), which is based 
on the most recent models for the nucleon-nucleon inter- 
action with the inclusion of a parameterized three-body 
force and relativistic boost corrections. Actually we adopt 
a simple analytic parameterization of this model given by 
[Heiselberg fc Hjorth- Jensen 1999| (later on HHJ). The lat- 
ter uses the compressional part with the compressibility 
K ~ 240 MeV, and a symmetry energy fitted to the data 
around nuclear saturation density, and smoothly incor- 
porates causality at high densities. The density depen- 
dence of the symmetry energy is very important since 
it determines the value of the threshold density for the 
DU process. The HHJ EoS fits the symmetry energy to 



Fig. 121 (left panel) shows the mass-radius relation 
and (right panel) the mass-central density relation. Solid 
lines are for HHJ EoS. In order to check an alterna- 
tive possibility and to demonstrate the sensitivity of 
the value of the DU threshold density to the selected 
model we also use a version of the relativistic non- 
linear Walecka (NLW) model in the parameterization of 
[Kolomeitsev fc Vos krescns ky 2003| The parameters of the 
NLW model are adjusted to the following bulk parame- 
ters of the nuclear matter at saturation: uq — 0.16 fm^'^, 
binding energy i?bind — —15.8 MeV, compression modu- 
lus K =250 MeV, symmetry energy agym = 28 MeV, and 
the effective nucleon mass m^(no) = 0.8 m^r. The corre- 
sponding coupling constants are as follows: 



ml 
91n t^n 



91n 



= 91.2506 
= 195.597 
= 77.4993 



b = 



0.00867497, 
0.00805981 , 



(6) 



(note that we treated obvious misprints in b and c in the 
paper of [Kolomeitsev fc Voskresensky 2003| ). This model 
was constructed in such a way that it reproduces almost 
the same thermodynamic properties as those in the model 
of |Heiselberg fc Hjorth-Jensen 1999| However, having not 
enough free parameters all NLW-based models yield signif- 
icantly lower threshold densities for the DU process than 
those given by variational calculations. [Migdal et al. 1990| 
treated this fact as a fragile point of the relativistic mean 
field models. Thereby they disregarded DU process from 
their subsequent analysis concentrating on the pion soften- 
ing and the PU possibility. ILattimer et al. 19911 used this 
fact to develop the DU-based scenario for the NS cooling. 
Anyhow, in the given NLW model the threshold density 
for the DU process is nj?" ~ 2.62 no (M°u ~ 1.29 Mq). 
Dots in Fig. (21 (right panel) indicate threshold densities 
for the DU process. An influence of the pion condensation 
on the EoS for n > n™ is assumed to be small and sup- 
pressed thereby. For the HHJ EoS the threshold density 



3 no corresponds to the NS mass M^^ ~ 1.32 M^ 



From Fig. 121 one can see that deviations in the M{n) re- 
lation for HHJ and NLW EoS are minor, whereas the DU 
thresholds are quite distinct. 

Fig. 121 demonstrates the concentrations of p, e and /Lt~ 
in HHJ and NLW models as a function of the baryon den- 
sity. All these dependencies are quite different for these 
two models. This means that relativistic mean field mod- 
els may describe well thermodynamic properties but yield 
quite different cooling picture compared with that given 
by more microscopically based variational calculations of 
Argonne-Urbana. The possibility of charged pion conden- 
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R [km] n(0) [fm'"^] 



the temperature profile being inhomogeneous during first 
10^ ^ lO'^ yr. Tlie influence of the crust on the coohng and 
heat transport is rather minor basically due to its rather 
low mass content. Thereby the temperature also changes 
slightly in the crust up to the envelope. 

3.3. Envelope 

Further on we need the relation between the crust 
and the surface temperature for NS. The sharp change 
of the temperature occurs in the envelope. This Tg — 
Tin relation has been calculated in several works, see 
IGlen fc Sutherland 1980l lYakovlev et al. 2003bl depend- 
ing on the assumed value of the magnetic field at the 
surface and some uncertainties in our knowledge of the 
structure of the envelope. Fig. 0] shows a range of avail- 



Fig. 2. Gravitational mass-radius relation and mass - 
central density dependence for the NS configurations cor- 
responding to the HHJ (solid lines) and NLW (dashed 
lines) model EoS. Dots indicate the DU threshold. 
Possibility of pion condensation is suppressed. 
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Fig. 4. Relation between the inner temperature and the 
surface temperature for different models used in our cal- 
culations. Dash-dotted curves indicate boundaries of pos- 
sible values = /(Tin), see lYakovlev et al. 2003bl 



Fig. 3. Densities of the charged particles as a function of 
baryon density for HHJ and NLW model EoS. Threshold 
density for the DU process is indicated. An influence of 
pion condensation for n > is neglected. 

sation is suppressed. Otherwise for n > the isotopic 
composition may change, see Migdal et al. 19901 in favor 
of an increase of a proton fraction and a smaller critical 
density for the DU reaction. 

3.2. NS crust 

The density n ~ 0.5 ^ 0.7 uq is the boundary of the NS 
interior and the inner crust. The latter is constructed of 
a pasta phase discussed bv lRav enhall et al. 19831 see also 
recent work of [Maruyama et al. 20041 Then there is the 
outer crust and the envelope. Note that our code generates 



able Tg — /(?!„) curves, taken from lYakovlev et al. 2003bl 
The solid curve, Tg — T^^ = /(Tin), corresponds to our 
best fit of the "slow cooling" objects. It matches the up- 
per boundary (the dash-dotted curve 77 — 4.0 • 10""'^^) for 
rather low Tg (in interval logTs[i<r] ~ 5.8 ^ 6) and the 
lower boundary (the dash-dotted curve = 4.0 ■ 10"^) for 
high Tg {\ogTs[K] > 6.8). The dash curve shows the so 
called simplified "Tsuruta law" T^'"''' = (10 Ti^f/^ used 
in many old cooling calculations. We will further vary dif- 
ferent possibilities. The dependence of the results on the 
Ts — Tin relation is demonstrated below in Figs.l^and llOl 
Figs. 112 d and d and Figs. UHl and CHI 

4. Main cooling regulators 

We compute the NS thermal evolution adopting our fully 
general relativistic evolutionary code. This code was orig- 
inally constructed for the description of hybrid stars by 
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IBlaschke et al. 20011 The main cooling regulators are the 
thermal conductivity, the heat capacity and the emissivity. 
In order to better compare our results with results of other 
groups we try to be as close as possible to their inputs for 
the quantities which we did not calculate ourselves. Then 
we add inevitable changes, improving EoS and including 
medium effects. 

4.1. Thermal conductivity 

We take the electron-electron contribution to the ther- 
mal conductivity and the electron-proton contribution for 
normal protons from lGnedin fc Yakovlev 19951 The total 
contribution related to electrons is then given by 

l/^e = 1/Kee + l/^ep- (7) 

For T > Tcp (normal "n" matter). For 
T < Tcp (superfluid "s" matter V lG^edin fc Yakovlev 19951 
suggested to drop the superfluid contribution 1/^1^. We 
use the expression 

^ep — ^ep/£,p > l^ep' i^) 

that gives a crossover from the non-superfluid case to the 
superfluid case. The vanishing of K|p for T <C Tcp is a 
consequence of the scattering of superfluid protons on the 
electron impurities, see lBlaschke et al. 20011 Following 
we get k" < k|. Surprisingly, it is in disagreement with 
Fig. 4 of IBaiko e t al. 2 0011 where probably the curves 
"SFl" and "SFO" should be interchanged. 
For the nucleon contribution, 

Kn = 1/ Knn + 1/ l^np j (9) 

we use the result of IBaiko et al. 20011 that includes cor- 
rections due to the superfluidity. Although some medium 
effects are incorporated in this work, the nucleon- hole cor- 
rections of correlation terms and the modification of the 
tensor force are not included. This should modify the re- 
sult. However, since we did not calculate k„ ourselves, we 
may only roughly estimate the modification. As we have 
shown above in Fig. 1, not too close to the critical point of 
the pion condensation the squared matrix element of the 
NN interaction pp^T^/uj'^ is of the order of 

the corresponding quantity |M|^j^j. ^ p"^ jj[ml. +p|-„] es- 
timated with the free one pion exchange, whereas iMj^^^ 
may significantly increase for n . To simulate the 

effect we just allow for the variation of k„ multiplying 
it by the factor = 10 (fast transport) and Ck — 0.3 
(slow transport), see Figs. 171 ITTl below. The former case 
is relevant for rather massive stars, whereas the latter, for 
rather light stars. A suppression of the nucleon-nucleon 
amplitude compared to the one for the free one pion ex- 
change model for n < ng is motivated by the in-medium 
T matrix calculations bv IBlaschke et al. 19951 

The total thermal conductivity is the straight sum of 
the partial contributions 

Ktot = Ke + Kn + ... (10) 

Other contributions to this sum are smaller than those 
presented explicitly [k^ and k„). 



4.2. Heat capacity 

The heat capacity contains nucleon, electron, photon, 
phonon, and other contributions. The main in-medium 
modification of the nucleon heat capacity is due to the 
density dependence of the effective nucleon mass. We use 
the same expressions as lSchaab et al. 19971 The main reg- 
ulators are the nucleon and the electron contributions. For 
the nucleons (i — n,p), the specific heat is l|Maxwell 1979|l 

c, - 1020(m*/m,) {nJnoy/\u Tg erg cm-^K"! , (11) 



for the electrons it is 
Ce - 6 X 10^9 (ne/no)'/' T 



9 erg cm 



K" 



(12) 



Near the phase transition point the heat capacity ac- 
quires a fiuctuation contribution. For the first order 
pion condensation phase transition this additional con- 
tribution contains no singularity, in difference with 
what would be for the second order phase transition, 
see IVoskresensky fc Mishustin 1982| |Migdal et al. 1990| 
Finally, the nucleon contribution to the heat capacity may 
increase up to several times in the vicinity of the pion 
condensation point. The effect of this correction on global 
cooling properties is rather unimportant. 

The symmetry of the 3P2 superfiuid phase allows for 
the contribution of Goldstone bosons (phonons): 



6 • lO'^T^ 



erg 
cm-^ K' 



(13) 



for T < Tcn{3P2), n > nc„(3P2). We also include this 
contribution in our study although its effect on the cooling 
is rather minor. 



4.3. Emissivity 

We adopt the same set of partial emissivities as in the 
work of I Schaab et al. 19971 The phonon contribution to 
the emissivity of the 3P2 superfluid phase is negligible. 
The main emissivity regulators are the MMU, see above 
rough estimation jSl, nPBF and pPBF processes, see 
above rough estimation 

Only qualitative behavior of the interaction shown in 
Fig. n is motivated by microscopic analysis whereas ac- 
tual numerical values of the correlation parameter and the 
pion gap are rather uncertain. Thereby we vary the values 
T{n) and w*^(n) in accordance with our above discussion 
of Fig. By that we check the relevance of alternative 
possibilities: a) no pion condensation and a saturation of 
the pion softening with increasing density, b) presence of 
pion condensation. 

We also add the contribution of the DU for n > n^^ , 
see above rough estimation 

All emissivities are corrected by correlation effects. 
The PU process contains an extra factor compared 
to the DU process. Another suppression of PU emissivity 
comes from the fact that it is proportional to the squared 
pion condensate mean field Near the critical point 

iy9p ~ 0.1 increasing with the density up to |(^|^ ~ /^/2, 
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where /^r ~ 93 MeV is the pion decay constant. Finally, 
the PU cmissivity is about 1-2 orders of magnitude sup- 
pressed compared to the DU one. Moreover, we adopt the 
same gap dependence for the PU process as that for the 
DU process. 




Fig. 5. Neutron and proton pairing gaps according to 
lYakovlev at al. 2003a] (thick solid, dashed and dotted 
lines) and to [Takatsuka fc Tamagaki 2004| (thin lines). 



4.4. Nucleon superfluidity 

In spite of many calculations, values of nucleon gaps in 
dense NS matter are poorly known. This is the conse- 
quence of the exponential dependence of gaps on the 
poorly known NN interaction for n ^ uq. Recent calcu- 
lations of lSchwenk fc Friman 20031 who included medium 
effects in evaluation of the 3P2 gap demonstrate its strong 
suppression up to values ^ 10 keV. Together with previous 
findings of ISchulze et al. 19961 ILombardo fc Schulze 20'00l 
they motivate to consider possibility of rather suppressed 
gaps. The suppression of the 3P2 gap for n > n^^ 
(for charged tt condensation), as explained in text of the 
work of lSchaab et al. 19971 is also in coincidence with the 
argumentation of Takatsuka fc Tamagaki 1997[ see also 
ITsuruta et al. 2002L ITsuruta 20041 

Below we start with the model used by 
lYakovlev at al. 2003al Also we use the gaps from 
recent work [Takatsuka fc Tamagaki 2004| to compare 
the results. The gaps are presented in Fig. |S1 Thick 
lines, the gaps from lYakovlev at al. 2003ai and thin 
lin es, from [Takatsuka fc Tamagaki 2004| the model AV18 
by [Wiringa et al. 1995[ The ISq neutron gap is taken 
the same in both models. We allow for a variation of 
gaps in wide limits in order to check the sensitivity of 
the results to their values. As we will see the cooling 
curves are sensitive to the magnitudes and to the density 
dependence of the gaps. Therefore further microscopic 
studies of the gaps are required. 



4.5. KDU, HDU, QDU and other 

The phase structure of dense NS matter might 
be very rich, including n^, tt^ condensates and 
K'^, K~ condensates in both S and P waves 
| [Kolomeitsev fc Voskres ensky 2003') ; charged p-meson 
condensation by fVoskresensky 1997[ ) ; couphng of 
condensates [Umeda et'^T7T994| fermion condensation 
yielding an efficient DU-like process in the vicinity 
of the pion condensation point (with the cmissivity 
- 10^7 T|, TO^ ( x l/T), see [Voskresensky et al. 2000| 
hyperonization, see [Takatsuka fc Tainag akr2004': quark 
matter with different phases, like so called 2SC, CFL, 
CSL, plus their interaction with meson condensates, see 
[K. Rajagopal fc F. Wilczek 2000[ IBlaschke et al. 20011 
and refs therein; and different mixed phases. In the 
present work, we suppress all these possibilities of extra 
efficient cooling channels. They are effectively simulated 
by our PU choice. The quark matter effects need a special 
discussion. We will return to the latter possibility in a 
subsequent publication. 

5. Numerical results 

5.1. Cooling of normal NS 

Although we have no doubts about the presence of nu- 
cleon superfluidity in NS interiors, we consider first the 
case of the complete absence of nucleon superfluidity. The 
reasons for that are as follows: i) Thus we compare our 
results with previous calculations, ii) On this example we 
select more essential and less essential ingredients in this 
many-parameter problem. Results are more transparent 
being demonstrated on a simplified example since the gen- 
eral case introduces new uncertain parameters and is more 
involved, iii) The actual values of gaps might be essen- 
tially smaller than those estimated in the literature, cf. 
ILomhardo fc Schulze 2000llS^hwenk fc Friman 20031 This 
is because most calculations did not include a proper 
medium dependence of the NN interactions. Thus we dis- 
cuss the limiting case of largely suppressed gaps. 

Fig. demonstrates the cooling evolution (for T^*'{t) 
dependence) of a normal NS for the HHJ EoS. Medium 
effects and PU possibility are disregarded ( "Standard sce- 
nario"). We see that one could easily explain the "slow 
cooling" points. Curves for NS masses in the range M ~ 
(1 1.839) Mq he in this region. For M ~ 1.839M0 the 
efficient DU process is switched on within the HHJ EoS 
{Standard + exotics DU" scenario). Thereby, the curve 
corresponding to M ~ 1.841Mq jumps down and already 
explains the "rapid cooling" points. The "intermediate 
cooling" points can be explained either by a very low NS 
mass (see curve for M ~ 0.5 Mq, or by M ~ 1.840 Mq. 
However, it seems rather unrealistic that 4 from 10 objects 
either relate to very low masses like M ~ 0.5 Mq, or are 
highly massive, as M ~ 1.840 Mq, where the DU process 
is switched on. "Rapid cooling" points are explained by 
very massive objects {M ~ 1.841 Mq) with DU process. 
They could be also explained by very low-mass objects 
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Fig. 6. Cooling of HHJ configurations without medium 
effects and PU for different masses of NS, Tg — Ti„ relation 
according to our fit. 

(M ~ 0.1 Mq). We further drop the latter possibility as 
rather unrealistic one and consider M > 0.5 Mq. 

Thus the picture as a whole looks unsatisfactory. Only 
"slow cooling" data are appropriately explained in a rea- 
sonable NS mass interval 1.0-^-1.839 Mq. The explanation 
of "intermediate cooling" and "rapid cooling", although 
possible, needs very unnatural assumptions. 
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Fig. 7. The influence of the heat conductivity on the sce- 
nario of Fig. 13 Two representative configurations with 
masses M = 1.40 Mq (sohd hues) and M = 1.841 Mq 
(dashed lines). We vary the heat conductivity by factors 
C = 10 and C = 0.3. 

Fig[7|demonstrates the sensitivity of the cooling curves 
to the variation of k„ . We scale k„ by a factor ^ = 10 and 
by Q = 0.3. The former case is meaningful for heavy ob- 
jects whereas the latter one for low-mass objects. We see 
that both increasing and decreasing of the thermal conduc- 



Fig. 8. Cooling evolution of the NS with normal matter 
HHJ EoS (for T^*^) including the medium modification of 
MU and other processes (MMU, MNB, etc.), without pion 
condensation. 




Fig. 9. Same as Fig. |H1 including pion condensation for 
rt > 3 710. 

tivity does not change the picture as the whole, as well as 
the conclusion drawn above. Transport is relevant only up 
to the first 10"^ y when the details of the Tg — rela- 
tion as well as temperature and density dependences of 
the cooling regulators determine the response of the cool- 
ing curves to a rescaling of the heat conductivity by the 
factor 

Fig. |H1 allows for medium effects the strength of which 
is assumed to be saturated with increasing density (no 
TT condensation for n > 3 no). Thus we use curves la, 
lb of Fig 1. We see that medium effects being included 
in calculation of the emissivity (MMU, MNB and oth- 
ers) significantly decrease all the curves. This allows us to 
explain reasonably well "rapid cooling" points but "slow 
cooling" and "intermediate cooling" cannot be addressed 
without inclusion of superfluidity. A partial suppression of 
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Fig. 10. Same as Fig. |51 but with Tg — Tin relation given 
by the Tsuruta law. 

medium effects by scaling of lo* does not allow to improve 
the picture. 

Comparison of Fig.|Sland Fig. shows that inclusion 
of medium effects regulates the mass dependence of the 
curves. In Fig. the curves rise with increase of the NS 
mass (for objects with M < 1.89 Mq, below DU thresh- 
old) whereas in Fig.|Slthe trend is changed to opposite. 

Fig. O additionally allows for tt condensation for n > 
3no (M > 1.32 Mq for HHJ EoS). Now we use curves 
la, 2, 3 of Fig 1. The picture remains almost the same as 
that shown in Fig.|Hl However objects with M > 1.32 Mq 
are cooled still faster by the efhcient PU reaction. If we 
assumed n^^ ~ 2.5 uq the pion condensation would start 
for NS masses M > 1.08 Mq. 

Fig.llOldemonstrates the same cooling evolution of nor- 
mal NS as Fig. ini but for the Tsuruta relation T/'^'Xi)- 
One can explain the "rapid cooling" and the "interme- 
diate cooling" but cannot explain the "slow cooling". We 
checked that the use of other Tg— Tin relations in the frame 
given by Fig. 01 does not change general trends. Medium 
effects and the possibility of tt condensation, although they 
regulate the NS mass- dependence show too rapid cooling. 
We see that without superfluidity the picture is unsatis- 
factory. Thereby we conclude that the cooling data call for 
the nucleon superfluidity. 

5.2. Cooling of superfluid NS 

First we check the best fit of the model 
IKaminker et al. 20011 lYakovlev at al. 2003al where 
gaps are given by Fig. El (thick lines) . In order to get a 
fit of the data within their "Standard + DU" scenario 
IKaminker et al. 200 11 lYakovlev at al. 2003al were forced 
to additionally switch off the 3P2 neutron gap. Since 
the full switching off the 3P2 gap is not supported by 
microscopic calculations we simulate the same effect by 
introducing the scaling factor 0.1 for the 3P2 neutron 



gap shown in Fig. |S1 by corresponding solid line. Then 
the magnitude of the 3P2 gap becomes to be com- 
parable with the value following from the findings of 
ISchwenk & Fri man 20031 who estimated the 3P2 neutron 
gap including medium effects. 
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Fig. 11. Cooling of NS configurations with superfluid nu- 
clear matter without medium effects and pion condensa- 
tion. The gaps are taken as in lYakovlev at al. 2003al see 
Fig. |31 the neutron 3P2 gap is additionally suppressed by 
a factor 10, - Tin relation is given by Tf*(t). 

Fig. II II demonstrates the cooling evolution {T^*^{t) de- 
pendence) of the above choice of gaps for a NS with 
HHJ EoS. Medium effects are not included. The curves 
rise compared to those in Fig. 6. We see that for M < 
1.839 Mq when the DU process is switched off, all the 
curves demonstrate very slow cooling even not explain- 
ing the "slow cooling" points (for M > 1 Mq). We also 
checked the artificially suppressed rate of the pPBF pro- 
cess used by IKaminker et al. 20011 lYakovlev at al. 2003al 
Suppression of the rate by a factor of 10 does not sig- 
nificantly affect the curves and does not change the con- 
clusion. The general trends of the curves are similar to 
those of IKaminker et al. 20011 However within the HHJ 
EoS the "intermediate cooling" points and "rapid cool- 
ing" points can be explained only by objects with M > 
1.839 Mq and DU cooling what seems us quite unsat- 
isfactory. Comparison of Fig. ^] with Fig. (BJ computed 
within the same scenario but without inclusion of the su- 
perfluidity, shows that both choices suffer of the very same 
shortcomings. The data call for medium effects. 

Fig.ElaUows for the medium effects (we use curves la, 
lb of FigQ] no tt condensation for n > 3 no). We smoothly 
cover all data points. All the data are explained by masses 
in the interval M ~ 1.39 1.84M0. The regular behavior 
is clearly seen. More massive objects cool faster than less 
massive ones. The star with the mass 1.25 Mq may relate 
to the old and hot objects like PSR 1055-52. 
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Fig. 12. Same as in Fig. 1111 including medium effects Fig. 14. Same as in Fig.lT^but T„ — Ti„ relation according 
without pion condensation (curves la, lb in Fig. to lYakovlev at al. 2003al r] = 4 x 10~^^. 





iog(t[yri) 



Fig. 13. Same as in Fig. El but using the Tsuruta law. Fig. 15. Same as in Fig. IT^ including tt condensation. 



We checked a sensitivity of the result to the depen- 
dence Ts — Tin. Fig. Pluses the "Tsuruta law" whereas 
Fig- El probes the choice of lYakovlev at al. 2003al rj = 
4.0 X 10^^^. We see that a variation of Tg — Tin does 
not change the picture as the whole. Only the interval 
of masses that cover the data is slightly changed. This in- 
terval is M ~ 1.36 ~ I.SOMq according to Fig. [T^ and 
M ~ 1 ^ 1.75 Mq according to Fig.lTl 

The mass interval that covers the data is M ~ 1.37 
1.46 Mq. Fig. 1151 allows for tt condensation for n > 3 no. 
Now we use curves la, 2, 3 of Fig^ The Tg — Tin relation 
is according to our fit. 

Fig. 1161 shows the same as Fig. llSl but for the Tsuruta 
law. In both cases the picture remains almost the same 
as that shown in Fig. El and Fig. El respectively. The 
difference starts for masses > 1.32 Mq due to switching 
on the efficient PU process. The interval of masses that 
cover the data is M ~ 1.36 1.452 Mq. Using of other 



Ts — Tin dependences does not change the picture. Thus 
we may conclude that the tt condensation does not con- 
tradict the data but, on the other hand, this assumption 
is not motivated by the existing cooling data. The hy- 
peron enhanced cooling or kaon condensation cooling do 
not change the picture if their critical densities are ^ 3 tiq. 
Note that the transition density can't be too low. At least 
n^^ > 2.5 ~ 2.7 no for the charged pion condensation 
within given EoS. The neutral PU process might be ad- 
ditionally suppressed compared to the charged one, see 
[Voskresensky 2001] that may weaken the above restric- 
tion. Otherwise we would get too fast cooling already for 
low mass objects and the NS with M ~ 1.4 Mq would 
cool faster than it is available by modern data. 

Fig. 1171 demonstrates the sensitivity of the curves pre- 
sented in Fig.Elto the variation of k„. We scale k„ by the 
factor C = 10 and C = 0.3. We see that both increasing 
and decreasing of the thermal conductivity k„ does not 
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Fig. 16. Same as in Fig.^lusing the Tsuruta law. 




Fig. 17. The influence of a change of the heat conductiv- 
ity on the scenario of Fig. El 

change the picture as the whole as well as the conclusion 
drawn above. As in the case when medium effects are sup- 
pressed, see Fig. [7| the effect of the transport is relevant 
for first 10^ yr. 

Figs. El - El show the dependence of the results on 
different variations of the gaps, their absolute values and 
the density dependencies. 

We remind the reader that the proton gap is artificially 
enhanced by 'Yakovl ev at al. 2003al Fig. El presents the 
same as Fig. 1151 but now we suppress the proton gap by 
a factor 0.5. An appropriate fit is achieved. In Fig. El 
we further suppress the proton gap (now by factor 0.2) 
and we also suppress the ISo neutron gap by factor 0.5 
simulating the medium effects in gaps. We again obtain an 
appropriate overall fit of the data. The NS masses covering 
available data are 1.33 1.44 Mq in case of Fig. 1181 and 
1.23-;-1.42 Mq in case of Fig. El However we have checked 
that scaling of all gaps shown by thick lines in Fig. Oby 
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Fig. 18. Same as Fig.Elbut with suppressed proton gap 
by a factor 0.5. 
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Fig. 19. Same as Fig.Elbut with 15*0 neutron gap sup- 
pressed by a factor 0.5 and 15*0 proton gap by a factor 0.2. 
The neutron 3P2 gap is remained to be suppressed by 0.1. 



the factor 0.1 already does not allow for the appropriate fit 
of the data. Thus we demonstrated important role played 
by all three types of the pairing: ISq nn and pp and 3P2 
nn. 

In Fig. El we include the same medium effects as 
those in Fig. El but we use the gaps according to 
[TakatsuEa fc Tamagaki 2004[ see thin lines in Fig. [3 ad- 
ditionally suppressing 3P2 gap by a factor 10. We see that 
with HHJ EoS and the given choice of the gaps we ob- 
tain the best fit of the whole set of the data. Thus the 
pion softening and appropriate density dependence of the 
gaps are quite sufficient to explain the modern NS cooling 
data. The data are covered by NS in wide mass interval 
0.5 1.84 Mq. For masses 1.0 1.84 Mq the picture is 
quite regular. Less massive stars cool slower, more massive 
stars cool faster. 
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Fig. 20. Same as in Fig. E| using the gaps of 
[Takatsuka fc Tamagaki 2004| see Fig. [SJ with additional 
suppression of the 3P2 gap by a factor 10. 



Fig. 22. Same as Fig. II 51 without suppression of 3P2 neu- 
tron pairing gap. 
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Fig. 21. Same as in Fig.lJOlbut including pion condensa- 
tion. 



Fig. 12] presents the same as Fig.|5nibut now including 
the pion condensation. The picture as the whole remains 
the same. Only cooling of massive stars is touched that 
slightly narrows the appropriate NS mass interval from 
1.0 ~ 1.75 M© to 1.0 1.70 Mq. 

Fig. El and Fig. US present the same as Fig. 1151 and 
Fig. 1^ but now we take off the suppression factor 0.1 
for the 3P2 neutron gap. From Fig.|23we see that 2 'slow 
cooling" and 2 "intermediate cooling" points related to old 
NS are not explained. In case of Fig. [221 only two lowest 
data points are explained. We again support the statement 
that the suppression of the SP2 gap is indeed called for by 
the cooling data. 



Fig. 23. Same as Fig. 1211 without suppression of 3P2 neu- 
tron pairing gap. 

6. Conclusion 

We have shown that the most up-to-date observed NS 
cooling data can well be explained. Actually we deal with 
a many-parametric problem that allows for a variation of 
many quantities. Besides the data points might be par- 
tially shifted from their positions shown in figures due 
to a number of uncertainties and assumptions used at 
their analysis. By above figures we have illustrated dif- 
ferent possibilities discriminating more probable expla- 
nations from less probable ones. We elaborated the ex- 
ample of the probably most realistic and microscopically 
supported EoS of the V18 + Sv + UIX*. Actually we 
used a simple parameterization of this EoS suggested 
by [Heiselberg fc Hjorth- Jensen 1999| (HHJ). In this EoS 
the DU process does not show up to M ~ 1.839M0. 
Thus within this model the "Standard + DU" scenario 
would demonstrate that the majority of experimentally 
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measured cooling points relate to very massive NS, with 
M > 1.84 Mq. From our point of view such a scenario 
seems unrealistic. 

We exploit in-medium effects in the calculation of the 
emissivities and the pairing gaps, and, that is less impor- 
tant, in specific heat and heat conductivity. In general, 
medium effects result in a significant suppression of the 
superfluid gaps, especially of the 3P2 nn gap, and they 
enhance the cooling rates of the MMU and the MNB pro- 
cesses through the pion softening effect with the increase 
of the density. Without the latter effect the "Standard 
+ PU" scenario suffers from an internal inconsistency. 
Pion condensation cannot take place without prelimi- 
nary softening of the pion mode at lower densities. And 
V.S., recent argumentation for the pion condensation, cf. 
lAkmal et al. 19981 ISuzuki et al. 19 99' further motivates 
the presence of a precursor pion softening, details see in 
[Migdal et al. 1990| Pion condensation at n > 3no does 
not contradict the data but the data can be explained 
also without pion condensation and other so called "ex- 
otics" (KU, HDU, etc) but with the pion softening. Based 
on the works ISchaab et al. 19971 IBlaschke et al. 20011 we 
have further improved our code. We showed that the re- 
sults are sensitive to the values and the density dependen- 
cies of the nucleon superfluid gaps. The strong suppression 
of the 3P2 go-P motivated by theoretical evaluations that 
incorporate medium effects is indeed required for the fit 
of the data. We show that all three groups of points "slow 
cooling", "intermediate cooling" and "rapid cooling" are 
now well explained on the basis of the "Nuclear medium 
cooling scenario" demonstrated here, where an important 
role is played by in-medium effects, cf. Voskresensky 2001, 

We may draw the following main conclusions. 

i) The normal matter assumption (see Figs. - I10|l 
seems rather unrealistic, as by itself, as in relation to 
the cooling data. One could explain the data but at the 
price of ignoring of medium effects in MMU and MNB. 
Then the "intermediate cooling" points can only be ex- 
plained by very low NS masses (as 0.5 Mq) or, together 
with "rapid cooling" points, by the very high NS masses 
M > 1.839 Mq, which allow for the DU process. In the 
latter case the mass window separating the "intermediate 
cooling" and "rapid cooling" is very narrow. Above price 
seems us too high and we drop such a scenario. The su- 
perfluidity is called for by the data. 

ii) Including superfluid gaps we see, in agreement with 
recent microscopic findings, that 3P2 neutron gap should 
be as small as ^ 10 keV. Otherwise one cannot explain at 
least several old objects (see Figs. 1221 123f) . Proton and ISq 
neutron gaps also might be suppressed by factors of order 
of several, as it is demonstrated by Figs. ^J^J [5011^ but 
the suppression by factor > 10 is not already permitted. 

iii) Medium effects associated with the pion soften- 
ing are called for by the data. As the result of the pion 
softening the pion condensation may occur for n > n^^ 
(n > Shq in our model). Its appearance does not contra- 
dict to the data (see Fig. [T5|l but also the data are well 
described, if the softening effect is rather saturated (see 



Fig. I12|) with increase of the density (as demonstrated by 
curves la, lb in Fig.^. At the same time the critical den- 
sities for the efficient DU-like exotic processes (such as PU 
on charged pion) should not be as small as < 2.5 no, if the 
given HHJ EoS is indeed correct. Otherwise one would get 
too rapid cooling of the NS of the typical 1.4 Mq mass. 
This also means that the proper DU threshold density 
can't be too low that puts restrictions on the density de- 
pendence of the symmetry energy. Both statements might 
be important in the discussion of the heavy ion collision 
experiments. 

iv) We demonstrated a regular mass dependence: for 
the NS masses M > 1.0 M© less massive NS cool slower, 
more massive NS cool faster. 

As we have mentioned, for the sake of simplicity we did 
not include the possibility of the hyperonization and other 
possibilities, like kaon condensation and fermion conden- 
sation, which may stimulate a more rapid cooling, being 
working in a line with the pion condensation. We did not 
include possible quark effects. The latter need a special 
treatment. The possibility of the color superconductivity 
in dense NS interiors opens a number interesting possibil- 
ities like the so called two-flavor color superconductivity 
(2SC) phase, color-flavor-locking (CFL) phase, color-spin- 
locking (CSL) phase. Their cooling is essentially different. 
We will return to this discussion in the nearest future. 
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